Dependence of the dielectric constant of electrolyte solutions on ionic concentration 
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We study the dependence of the static dielectric constant of aqueous electrolyte solutions upon the 
concentration of salt in the solution and temperature. The model takes into account the orientation 
of the solvent dipoles due to the electric field created by ions, the ionic response to an applied field, 
and the effect of thermal fluctuations. The analysis suggests that the formation of ion pairs by a 
small fraction of disassociated ions can have a significant effect on the static dielectric constant. The 
model predicts the dielectric has the functional dependence e(c) = e w — f3L(3ac/ f3) where L is the 
Langevin function, c is the salt concentration, e w is the dielectric of the pure water, a is the total 
excess polarization of the ions and /3 is the relative difference between the water dipole moment and 
the effective dipole moment of ion pairs as weighted by the density of ion pairs and their structural 
rigidity. The functional form gives an extremely accurate prediction of the dielectric constant for a 
variety of salts and a wide range of concentrations by fitting only the parameter /3. 

PACS numbers: 61.20.Qg,77.22.Gm,77.22.Ch 



It is difficult to overstate the importance of aqueous 
electrolyte solutions in biological and electrochemical sys- 
tems. There has been extensive study of physical prop- 
erties of such solutions over the last 120 years. The first 
systematic experimental study of the dielectric properties 
of salt-water solutions was conducted in 1948 by Hasted 
et al. pQ. In this work, the static dielectric constant of a 
solution was observed to decrease with the salt concen- 
tration, a phenomena called dielectric decrement. Intu- 
itively, the dielectric decrement stems from the fact that 
the local electric field generated by each ion dominates 
the external applied field. The polar water molecules 
tend to align with the local ionic field, creating a hydra- 
tion shell around the ion, lowering the response of the 
water molecules to the external field and hence lowering 
the dielectric constant. In dilute solutions (typically salt 
concentrations less than 1.5M) the dielectric decrement 
is linear, 

e = e w -ac, (1) 

where s w is the dielectric constant of pure water, c is 
the salt concentration, and a is a phenomenological ion- 
specific parameter, known as the total excess polariza- 
tion of the ionic species. At higher salt concentrations 
significant deviations from linearity are observed and the 
dielectric decrement is observed to saturate [T]. 

Haggis et al. [5] modeled the observed linear dielectric 
decrement by considering the hydration shells as small 
spherical regions with a low dielectric constant, immersed 
in pure water medium with a high dielectric constant. 
The macroscopic dielectric constant of the solution was 
then computed by homogenization. The model was later 
refined by considering the variation of the local dielectric 
constant near the ions [3] and finite-size effects [J. The 
treatment of hydration shells as spheres, however, is jus- 
tifiable only for dilute solutions (typically less than 1M) 
for which the hydration shells do not overlap. 



Only very recently, Levy, Andelman and Orland [3] 
presented a model for the dielectric decrement in con- 
centrated solutions. Using a field-theory approach they 
calculated the average dielectric constant around each ion 
at the mean-field level, and accounted for ion-ion inter- 
actions and hydration shell overlap via a one-loop correc- 
tion. The resulting prediction of the dielectric constant 
affords a good fit to data from a large range of concentra- 
tions of different salts using a single fit parameter related 
to the effective size of the ions. 

In this letter, we develop a model for high salt concen- 
trations (in excess of 1.5M) in which the solvation shells 
of the ions strongly overlap. The model is not based upon 
the field around a single ion. Rather, it assumes the wa- 
ter dipoles are influenced by an aggregate of ions, i.e., 
the local field acting on a water dipole arises from the 
net fields created by a surrounding ionic configuration. 
A key element in the analysis is based on the well known 
fact that a small fraction of the ions form ion pairs [BJ. 
The analysis in this letter suggests that the propensity 
of even a small fraction of the ions to form pairs has a 
significant, even dominant, effect on the static dielectric 
constant of the solution. The key result is an analytic 
prediction for the solution dielectric constant as a func- 
tion of salt concentration c, 

£ (c)= £w -/3i(^c), (2) 

where L is the Langevin function 

L(v) = coth(w) , 

v 

e w is the dielectric of the pure solvent, a is the total 
excess polarization of the ions and j3 is the relative dif- 
ference between an effective ion-pair dipole moment and 
the water dipole moment. This functional form gives 
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an extremely accurate prediction of the static dielectric 
constant over a large range of concentrations of different 
salts with a single fit parameter (3 for each mixture. 
The static dielectric constant takes the functional form 



e(c) 



dP(c, E CK ) 



0E o 



(3) 



E cx =0 



where E cx is external field intensity, P is the orientational 
polarization of the medium, and is the contribution 
to the dielectric due to other effects such as molecular 
polarization. Assuming is independent of the ionic 
concentration c, the aim of this letter is to suggest a 
mechanism for the dependence of the orientational po- 
larization P upon c. 

The first stage is to determine the contribution of a 
single water dipole to the orientational (dipolar) polar- 
ization, by introducing an ionic field to the standard 
Langevin dipole analysis, see, e.g., [7J sect. 4.6] or [8l p. 
2f4] for details. Consider a water dipole p surrounded 
by point-like ions under an applied external field E cx . 
The electric field due to the ions, in the absence of water 
dipoles, is denoted by Ei on . The potential energy W of 
the dipole p is given by 



W = -p • (E i( 



E n 



The energy W is minimized when the dipole is aligned 
with the field Ei on + E ox , hence the dipole orientation is a 
trade-off between the tendency of the dipole to align with 
the field E ion +E cx and thermal fluctuations that disrupt 
this ordering. The contribution of a given dipole p to the 
polarization due to the external field E cx is Proj Eox (p — 
Po) where po is the average orientation of the dipole p 
in the absence of an external field, i.e., in the direction 

of Ej nn , 



Ei. 



Po 



IE;, 



The expectation of the contribution of the dipole p to 
the external field E cx is given by the Boltzmann average 



P local — 



Je ^Proj Ecx (p-p ) 



(4) 



where E- lon — |Ei on |, and E cx = |E CX |. This integral can 
be explicitly evaluated following [5], 



local 



P L 



pE cx 
kT 



1 - L 



pE k 



kT 



(5) 



where p = |p|, see supplementary material for details. 

Expression ^ quantifies the contribution of a single 
water dipole to the orientational polarization. The sec- 
ond stage of the derivation considers the aggregate con- 
tribution P of all water dipoles to the orientational po- 
larization. Assuming the water dipoles are uniformly dis- 
tributed in space, P is given by 



(6a) 



P{c)= / f(E ion ;c)P local (E ion )dE k 



where Pi oca i is given by ^ and f(Ei on ; c) is the proba- 
bility density function of the ionic field intensity in the 
absence of water, given by [TO] , 



f{E ion ;c) 



4 f 
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P 
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Kn = <**c. 

(6b) 

Here, E* on is the most probable ionic field intensity 
and a*(M _1 ) is a dimensional constant of proportion- 
ality between the normalized ionic field intensity and the 
ionic concentration. The derivation of f(E- lon ) in [TU] re- 
lates a* to the electric field screening length. 

The effect of water on the ionic field intensity is ne- 
glected in these calculations, i.e., the ionic field is con- 
sidered as if the ions were in vacuum. The justification 
for this simplification is that in high concentrations, the 
water/ion ratio does not allow efficient screening by wa- 
ter molecules. For example, in 1M concentration, the 
average distance between ions and their counter-ions can 
not exceed 1.2nm. At this length scale, at most 3 water 
molecules can reside along the line segment connecting 
the ionic centers, hence the screening of electric field by 
the water is limited. In contrast, in ImM solutions, ion 
separation distances can average 12nm, with hundreds of 
water molecules engaged in screening. 

A key feature of the model is that a*, which re- 
lates ionic concentrations to screening length, particu- 
larly ionic screening, is independent of ionic concentra- 
tions above 1.5-2M. At these concentrations inter-ionic 
distances relate weakly to the concentration, and the abil- 
ity of ions to redistribute to improve effective screening 
becomes limited due to finite size effects and thermal 
fluctuations. As a result, the effective screening length 
saturates. 

The Langevin function can be approximated by 



L(v) w 1 - 1/v 



(7) 



for v ^> vl = 0.2. Using ^ to simplify ^ in the inte- 
gral ^ yields, 



P(c)^pL' pEcA 2 



kT 



(8) 



Although the primary goal of this derivation is to obtain 
a prediction for P(c) in the high concentration regime, it 
is still helpful to obtain an approximation that is compat- 
ible with known results in the dilute regime. Relation 
implies that the orientational polarization is linear for low 
concentrations. Using ^ in (JsJ) to obtain an approxima- 
tion of P that is linear for a*c <C I yields 



P(c) w pL 



pEcx 
kT 



1 - L 



(9) 



The resulting expression ([9| well approximates the tail 
of P in ^ and in addition complies with the observed 
linear decrement of the orientational polarization at low 
concentrations where ^ is invalid, see Figure [T] 



3 




FIG. 1. (color online) Orientational polarization (JsJ for a* = 
2 (-) and the approximation |9| (--). 

Expression ([9| quantifies the aggregate contribution of 
all water dipoles to the orientational polarization P. As 
expected, P — > as c — > oo (or £ ion — > oo), since the wa- 
ter dipoles become 'immobilized' by the strong ionic field. 
Experimental data, however, suggests that the dielectric 
constant of many electrolyte solutions approaches a lim- 
iting value, in some instances as high as e — 20, as the 
ionic concentration increases, see, e.g., Figure [2p. This 
observation supports the idea that water dipoles are not 
the sole contributors to the orientational polarization, 
but that ionic pairs may also contribute to the orien- 
tational polarization. This line of thought is supported 
by the study of molten salts in the absence of a solvent, 
in which orientational polarization is significant and the 
static dielectric constant is relatively high, e.g., e = 20 
for Ethylammonium Nitrate 

To form an effective dipole, ions need to tightly paired. 
Molecular dynamics simulations and small angle X-Ray 
scattering experiments conducted in a wide range of 
concentrations, 0.01M < c < 20M, consistently show 
that the locations of the peaks of the anion-cation ra- 
dial distribution functions are independent of concentra- 
tion, see, e.g., [T2TU7| . These anion-cation radial dis- 
tribution functions suggest that for c > 0.01M, a por- 
tion of the anions and cations associate into ion-pairs 
with ion-specific distances of r = 2.5 — 3.5A° for contact 
pairs or r = 4.5 — 5.5/1° for looser pairs. A fundamental 
assumption of the model is that these ion pairs act as 
dipoles and consequently contribute to the orientational 
polarization of the medium. As a first approximation, 
the effect of the ion-pairs is modeled by a local effective 
dipole moment pf^ ctlvc of uniform magnitude which is 
aligned in the direction of the local ionic field 



effective 
Pion 



(x) 



Eio n (x) 

|E ion (x)| 



Pion /"ion -Sion 



poffecth 



Here, p lon is the dipole moment of an ion-pair, /; on is the 
fraction of ions that form ion-pairs and < s- lon < 1 is 
the structural rigidity of the dipole assumed to be high 
for contact pairs, and lower for looser pairs. 

Incorporating the contribution, Proj E pf^ ctlvc , of the 
ion-pair dipole to the orientation polarization into the 



Boltzmann averaging Q and evaluating the integral 
yields 



Pi 



local 



pL 



pE cy . 
kT 
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kT 



(10) 



where 



p-i'i,, 



plify ( 10 ) in the integral Q yields 
fpE c 



As before, using ^ to sim- 



P(c) = pL 



\ kT 



ira c 
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(11) 



Finally, substituting (11) into ([3| gives the functional 
form ([2]) where the empirical parameters in take the 
form 



(i 



(3a* (3 = (e v 



The functional relation ^ is validated against six dif- 
ferent sets of experimental data, where the parameter /3 
is fitted separately for type of electrolyte. Figure [2j\ 
presents the experimental data of Hasted et al. [1] for 
an NaCl solution at 21°. Prediction ^ agrees very well 
with the experimental data over the full range of < c < 
6M. Figure [2^3 presents the experimental data of Wei et 
al. [18] for LiCl, RbCl and CsCl solutions at 25°. The 
predicted dielectric ^ remains accurate for LiCl data at 
concentrations as high as 13M. Finally, Figure [2p shows 
data taken from (18j , which were complied and presented 
in [5], together with the prediction obtained using the 
field-theory approach [5] (solid black). The functional 
form ([2]) for the dielectric permits a resolution of the dif- 
fering high-concentration dielectrics of the KC1 and KF 
salts, implying that KC1 has a greater predilection to 
form ion-pairs than KF. 

The model naturally incorporates the dependence of 
the dielectric constant upon temperature T, in which 
case all parameters become a function of temperature. 
The validity of functional relation ([2| is validated against 
four different sets of experimental data, due to Buchncr 
et al. 19J, for NaCl at various temperatures, where the 
parameters a and (3 are fitted separately for each tem- 
perature, see Figure [3] 

In conclusion, we have derived a model which gives 
an extremely accurate prediction for the static dielectric 
constant of aqueous electrolyte solutions. The key pa- 
rameters arising from the model are the constant of pro- 
portionality between the most probable ionic field inten- 
sity and the ionic concentration a* and the parameter (3 
related to the relative difference between the water dipole 
moment and the effective ionic dipole moment /^,n CCtlvc . 
Assuming that only rigid contact ion-pairs contribute to 
the orientational polarization, the fitted values of (3 in 
Figure[2]imply that / ion = 0.024-0.073, i.e., 2.4%-7.3% 
of the ions associate into ion-pairs, excluding the KF salt 
for which data suggests a weak propensity for K + -F~ to 
form ion pairs in water. These values are compared with 
an approximate of 8% ion-pairs in molten salt Ethylam- 
monium Nitrate [TT]. The fitted values of (3 in Figure [3] 
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FIG. 2. (color online) Comparison of the predicted dielectric constant KSJ, with experimental data as function of ionic con- 
centration c for various salts. The theoretical prediction was calculated using /3 as a fitting parameter. A: Data for NaCl salt 
from Q], compared to (2} with /3 = 42.07. B: Data from [18], where fit for RbCl and CsCl salt (--) is obtained with f3 = 67.20 

and for LiCl (-) with (3 = 53.22. C: Figure 2(b) from [5] where fit for KF ( ) is obtained with (3 = 75 and with (3 = 60 for 

KC1 (. — ). Solid black curve is the prediction obtained using the field-theory approach [5]. 
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FIG. 3. (color online) Comparison of the predicted dielectric 
constant |2| , with experimental data from |19| as a function of 
ionic concentration c for NaCl at various temperatures. Data 
for T = 5° fit with a = 13.7 and (3 = 70.25 (- -), for T = 20° 
with a = 12 and (3 = 52.94 (-), for T = 25° with a = 11.5 
and f3 = 47.91 (.-) and for T = 35° with a = 10.7 and /3 = 
40.12 (x-). 



show that /3 decreases with temperature, which is at- 
tributed to a decrease in the structural rigidity of the 
ionic pairs as temperature increases. 

A key implication of this work is a reinterpretation of 
the excess polarization parameter as the constant of pro- 
portionality between the most probable ionic field inten- 
sity and the ionic concentration, weighted by /?. This 
interpretation is consistent with the data showing an 
anomalously high excess polarization of H + ions. The 
value a = 20 for HC1 [T], compared to a = 11 — 14 



for the various salt solutions of Figure [2j reflects the 
high entropic cost associated with the localization of pro- 
tons [201121], leading to higher value of both a* and /3. 
Indeed, this derealization inhibits effective pairing of the 
ions, and also inhibits the effectiveness of protons in the 
screening of counter ions, as compared to K, Li, and Na. 
This suggests that the constant of proportionality a* be- 
tween the ionic field intensity and the ionic concentration 
should be significantly higher in salts which disassociate 
into H + . Similarly, we observe that the dimensional con- 
stant of proportionality between the ionic field intensity 
and the ionic concentration kTa* /p increases with tem- 
perature, since thermal fluctuations reduce the effective- 
ness of electric field [2"2"| . 

Electrolyte solutions are highly complex mixtures, and 
the model neglects numerous structural elements: wa- 
ter/water interactions, water/ion-pair interactions, finite 
size effects, reaction field effects, influence of the hydro- 
gen bond network, and the decrease in water molarity as 
ionic concentration increases. While a systematic study 
of these effects is clearly important, it is plausible that 
they will contribute perturbatively. 

Finally, utilizing the explicit evaluation of the Boltz- 
mann average Q , the model is readily extendable to the 
case of a strong applied field. In addition, although the 
model is validated only with 1:1 electrolytes for which rel- 
evant experimental data is available, it is also expected to 
be valid for 2:2 electrolytes. In this case, the anticipated 
higher tendency to form ion-pairs, would be manifested 
in small values of j3*, leading to a high dielectric con- 
stant e — e w (l — /?*) at high concentrations. Further 
more, (3* may become negative |23j leading to a non- 
monotonic dependence of the dielectric constant upon 
concentration. 
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